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Abstract: We present a review of the chaotic hypothesis and discuss its applications to 
intermittency in statistical mechanics and fluid mechanics proposing a quantitative defi- 
nition. Entropy creation rate is interpreted in terms of certain intermittency phenomena. 
An attempt to a theory of the experiment of Ciliberto-Laroche on the fluctuation law is 
presented. 

§1. Introduction. 

A general theory of non equilibrium stationary phenomena extending classical thermo- 
dynamics to stationary non equilibria is, perhaps surprisingly, still a major open problem 
more than a century past the work of Boltzmann (and Maxwell, Gibbs,...) which made 
the breakthrough towards an understanding of properties of matter based on microscopic 
Newton's equations and the atomic model. 

In the last thirty years, or so, some progress appears to have been achieved since the 
recognition that non equilibrium statistical mechanics and stationary turbulence in fluids 
are closely related problems and, in a sense, in spite of the apparently very different nature 
of the equations describing them they are essentially the same. 

The unifying principle, originally proposed for turbulent motions by Ruelle, [Ru78], in 
the early 970's, has been extended to statistical mechanics and eventually called the 
chaotic hypothesis", [GC95]: 

Chaotic hypothesis: Asymptotic motions of a chaotic system, be it a multi particle system 
of microscopic particles or a turbulent macroscopic fluid, can be regarded as a transitive 
Anosov system for the purposes of computing time averages in stationary states. 

It may be useful to make a few comments on how this is supposed to be interpreted. 
The conclusions that we draw here from the chaotic hypothesis are summarized in §13 
which might be consulted at this point. For a review on the subject seen from a different 
perspective see [Ru99a] 

§2. Meaning of the chaotic hypothesis. 

Anosov systems are well understood dynamical systems: they play a paradigmatic role 
with respect to chaotic systems parallel to the one harmonic oscillators play with respect 
to orderly motions. They are so simple, and yet very chaotic, that their properties are 
likely to be the ones everybody develops in thinking about chaos, even without having 
any familiarity with Anosov systems which certainly are not (yet) part of the background 
of most contemporary physicists. 



Informally a map a; — > Sa; is a Anosov map if at every point x of the bounded phase space M one 
can set up a local system of coordinates with origin at x, continuously dependent on x and covariant 
under the action of S and such that in this comoving system of coordinates the point x appears as 
a hyperbolic fixed point for S. The corresponding continuous time motion, when the evolution is 
a: — > Sja:, * 6 -R, requires that the local system of coordinates contains the phase space velocity x as one 
of the coordinate axes and that the motion transversal to it sees a; as a hyperbolic fixed point: note that 
a motion in continuous time cannot possibly be hyperbolic in all directions and it has to be neutral in 
the direction of x because the velocity has to be bounded if M is bounded, while hyperbolicity would 
imply exponential growth as either t— *+cxDort— »— oo. Furthermore there should be no equilibrium 
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In general an Anosov system has asymptotic motions which approach one out of finitely 
many invariant closed sets C\, . . . ,Cq each of which contains a dense orbit, :one says that 
the systems {Cj,St) are "transitive". One of them, at least, must be an attractive set. 

To say that "the asymptotic motions form a transitive Anosov system" means that 

(1) each of the sets Cj which is attractive is a smooth surface in phase space and 

(2) only one of them is attractive: 

The last "transitivity" assumption is meant to exclude the trivial case in which there 
are more than one attractive sets and the system de facto consists of several independent 
systems. 

The smoothness of Cj is a strong assumption that means that one docs not regard 
possible lack of smoothness, i.e. fractality, as a really relevant property in systems with 
large number of degrees of freedom. In any event one could consider (if necessary) 
replacing "Anosov systems" with some slightly weaker property like "axiom A" systems 
which could permit more general asymptotic motions. Here we adhere strictly to the 
chaotic hypothesis in the stated original form, [GC95]. 

§3. Basic implications of the chaotic hypothesis and relation with the ergodic 
hypothesis. 

The chaotic hypothesis boldly extends to non equilibrium the ergodic hypothesis: applied 
to equilibrium systems, i.e. to systems described by Hamiltonian equations, it implies 
the latter, [Ga98]. This means that if a Hamiltonian system at a given energy is as- 
sumed to verify the chaotic hypothesis, i.e. to be a transitive Anosov system, then for 
all observables F {i.e. for all smooth functions F defined on phase space) 



where is the Liouville distribution on the constant energy surface M, and (3.1) holds 
for almost all points x G M, i.e. for x outside a set J\f of zero Liouville volume on M. 
Being very general one cannot expect that the chaotic hypothesis will solve any special 
problem typical of non equilibrium physics, like "proving" the Fourier's law of heat 
conduction, the Ohm's law of electric conduction or the K41 theory of homogeneous 
turbulence. 

Nevertheless, like the ergodic hypothesis in equilibrium, the chaotic hypothesis accom- 
plishes the remarkable task of giving us the "statistics" of motions. If M is the phase 
space, which we suppose a smooth bounded surface, and t Stx is the motion starting 
at a; e M, the time average: 



of the observable F exists for x outside a set of zero phase space volume and it is 
a;-independent, thus defining the probability distribution hsrb via (3.2). 

Note, in fact, that the probability distribution ^J^sRB defined by the l.h.s. of (3.2) is 
uniquely determined (provided it exists): it is usually called the statistics of the motion" 
or the ^^SRB distribution'' associated with the dynamics of the system. 

To appreciate the above property (existence and uniqueness of the statistics) the fol- 
lowing considerations seem appropriate. 



points and the periodic points should be dense in piiase space. When the system has one or more 
(the so called "hysteresis phenomenon") attracting sets which do not occupy the whole phase space the 
chaotic hypothesis can be interpreted as saying that each attracting set is a smooth surface on which 
the time evolution flow (or map) acts as an Anosov flow (map). 




(3.1) 




(3.2) 
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3.3 



An essential feature, and the main novelty, with respect to equilibrium systems is that 
non conservative forces may act on the system: this is in fact the very definition of "non 
equilibrium system" . 

Since non conservative forces perform work it is necessary that on the system act also 
other forces that take energy out of it, at least if we wish that the system reaches a 
stationary state, showing a well defined statistics. 

As a consequence any model of the system must contain, besides non conservative forces 
which keep it out of equilibrium by establishing "/j!oW on it (like a heat flow, a matter 
flow, ...), also dissipative forces preventing the energy to increase indefinitely and forcing 

the motion to visit only a finite region of phase space. 

The dissipation forces, also called "thermostatting forces" , will in general be such that 
the volume in phase space is no longer invariant under time evolution. Mathematically 
this means that the divergence — cr(.r) of the equations of motion will be not zero and 

its time average Jj^ (j{y)^sRB{dy) (j+ will be positive or zero as it cannot be negative 
("because phase space is supposed bounded": see [Ru96]). 

One calls a system "dissipative" if (T+ > and we expect this to be the case as soon as 
there are non conservative forces acting on it. 

We see that if a system is dissipative then its statistics ushb must be concentrated on 
a set of zero volume in M: this means that ij.srb cannot be very simple, and in fact it 
is somewhat hard to imagine it. 

If the acting forces depend on a parameter E, "strength of the non conservative forces" , 
and for E = the system is Hamiltonian we have a rather unexpected situation. At 
E = the chaotic hypothesis and the weaker ergodic hypothesis imply that the statistics 
MsflB is equal to the Liouville distribution hl; but if S ^ 0, no matter how small, 
itiuill not be possible to express ^isrb via some density psiy) iu the form ^i.sRB{dy) = 
PE{y)fiL{dy), because psrb attributes probability 1 to a set J\f with zero volume in phase 
space {i.e. Hl{J^) = 0). Nevertheless natura non facit saltus (no discontinuities appear 
in natural phenomena) so that sets that have probability 1 with respect to psRB may be 
all still dense in phase space, at least for E small. In fact this is a "structural stahilit'i}'' 
property for systems which verify the chaotic hypothesis (see [Ga96c]) 

The above observations show one of the main difficulties of non equilibrium physics: the 
unknown hsrb is intrinsically more complex than a function pE{y) and we cannot hope 
to proceed in the familiar way we might have perhaps expected from previous experiences: 
namely to just set up some differential equations for the unknown pE{y)- 

Hence it is important that the chaotic hypothesis not only guarantees us the existence 
of the statistics fisRB but also that it does so in a "constructive way" giving at the same 
time formal expressions for the distribution psRB which should possibly play the same 
role as the familiar formal expressions used in equilibrium statistical mechanics in writing 
expectations of observables with respect to the microcanonical distribution ^l- 

For completeness we write a popular expression for Hsrb- If 7 is a periodic orbit in 
phase space, x-y a point on 7, T(7) the period of 7 then 

/ F{y)i,sRB{dy) = lim 6^2^ijML /° '-^^ (3.3) 

J E,:T(,)<Te-^o '^(^-.)-T(7) 



This is simple in the sense that it does not require, to be formulated, an even slight 
understanding of any of the properties of Anosov or hyperbolic dynamical systems. But 
in many respects it is not a natural formula: as one can grasp from the fact that it 
is far from clear that in the equilibrium cases (3.3) is an alternative definition of the 
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microcanonical ensemble (i.e. of the Liouville distribution /ii), in spite of the fact that 
in this case a = and (3.3) becomes shghtly simpler. 

To prove (3.3) one first derives alternative and much more useful expressions for iisrb 
which, however, require a longer discussion to be formulated, see [Ga99a] , [Ga86c] : the 
original work is due to Sinai and in cases more general than Anosov systems, to Ruelle 
and Bowen. 

§4. What can one expect from the chaotic hypothesis? 

In equilibrium statistical mechanics we know the statistics of the motions, if the ergodic 
hypothesis is taken for granted. However this hardly solves the problems of equilibrium 
physics simply because evaluating the averages is a difficult task which is also model 
dependent. Nevertheless there are a few general consequences that can be drawn from 
the ergodic hypothesis: the simplest (and first) is embodied in the '^heat theorem" of 
Boltzmann. 

Imagine a system of A'' particles in a box of volume V subject to pair interactions and to 
external forces with potential energy Wy, due to the walls and providing the confinement 
of the particles to the box. Define 



where the averages are taken with respect to the Liouville distribution on the surface of 
energy U. 

Imagine varying the parameters on which the system depends {e.g. the energy U and 
the volume V) so that dU, dV are the corresponding variations of U, V, then 



expresses the heat theorem of Boltzmann. 

It is a consequence of the ergodicity assumption, but it is not equivalent to it as it only 
involves a relation between a few averages {U,p,V,T), see [Bo66], [Bo84], [Ga99a]. Not 
only it gives us a relation which is a very familiar property of macroscopic systems, but 
it also suggests us that even if the ergodic hypothesis is not strictly valid some of its 
consequences might, still, be regarded as correct. 

The proposal is to regard the chaotic hypothesis in the same way: it is possible to imagine 

that mathematically speaking the hypothesis is not strictly valid and that, nevertheless, 
it yields results which are physically correct for the few macroscopic observables in which 
one is really interested in. 

The ergodic hypothesis implies the heat theorem as a general ("somewhat trivial") 
mechanical identity valid for systems of N particles with N = 1, 2, . . . , 10^^, . . .. For 
small N it might perhaps be regarded as a curiosity: such it must have been considered 
by most readers of the key paper [Bo84] who were possibly misled by several examples 
with N = 1 given by Boltzmann in this and other previous papers. Like the example of 
the system consisting of one "averaged" Saturn ring, i.e. one homogeneous ring of mass 
rotating around Saturn with energy U, kinetic energy T and "volume" V (improbably 
identified with the strength of the gravitational attraction!). But for = 10^'^ it is no 
longer a curiosity and it is a fundamental law of thermodynamics in equilibrium: which, 
therefore, can be regarded on the same footing of a symmetry being a direct consequence 
of the structure of the equations of motion, [Ga99a] appendices to Ch.l and Ch.9. It 
reflects in macroscopic terms a simple microscopic assumption {i.e. Newton's equations 
for atomic motions, in this case). 



T =average kinetic energy 

U =total energy 

p =average of dyWy 



(4.1) 



{dU +pdV)/T = exact 



(4.2) 
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No new consequences of even remotely comparable; importance; arc known to follow from 
the chaotic hypothesis besides the fact that it implies the validity of the ergodic hypoth- 
esis itself (hence of all its consequences, first of them classical equilibrium statistical 

mechanics). 

Nevertheless the chaotic hypothesis does have some rather general consequences. We 
mention here the fluctuation theorem. Let a{x) be the phase space contraction rate and 
a+ be its SRB average {i.e. a+ = J a{x) HsRsidx)), let r > and define 

p{x)=T-^ ^^dx (4.3) 

J-t/1 0-+ 

and study the fluctuations of the observable p{x) in the stationary state hsrb- We write 
'Kr{p)dp the probability that, in the distribution Hsrb, the quantity p{x) has actually 
value between p and p + dp as 

iTr{p)dp = const e''^^^^'^dp (4.4) 
Them lim^r^oo Ct{p) = C(p) exists and is convex in p; and 

Theorem: (fluctuation theorem) Assume the chaotic hypothesis and suppose that the 
dynamics is reversible, i.e. that there is an isometry I of phase space such that 

ISt = S-J, = 1 (4.5) 

and that the attracting set is the full phase space.^ Then 

C{-p) = Cip)-^+P, for all p (4.6) 

where a+ = hsrb{(^)- 

It should be pointed out that the above relation was first discovered in an experiment, 
see [ECM93], where also some theoretical ideas were presented, correctly linking the 
result to the SRB distributions theory and to time reversal symmetry. Although such 
hints were not followed by what can be considered a proof, [CG99], still the discovery 
has plaid a major role and greatly stimulated further research. 

The interest of (4.6) is that, in general, it is a relation without free parameters. The 
above theorem, proved in [GC95] for discrete evolutions (maps) and in [Ge97] for con- 
tinuous time systems (flows), is one among the few general consequences of the chaotic 

hypothesis, see [Ga96a], [Ga96b], [Ga99b] for others. 

§5. Non equilibrium ensembles. Thermodynamic limits. Equivalence. 

The chaotic hypothesis gives us, unambiguously, the probability distribution fisRS which 
has to be employed to compute averages of observables in stationary states. 

For each value of the parameters on which the system depends we have, therefore a 
well defined probability distribution hsrb- Calling a = {ai, . . . ,ap) the parameters 
and iia the corresponding SRB distribution we consider the collection £ of probability 



It is perhaps important to stress that we distinguish between attracting set and attractor. the first is a 
closed set such that the motions that start close enough to it approach it ever closer; an attractor is a 

subset of an attracting set that 

(1) has probability 1 with respect to the statistics ju of the motions that are attracted by the attracting 
set (a notion which makes sense when such statistics exists, but for a zero volume set of initial data, 
and is unique) and that 

(2) has the smallest Hausdorff dimension among such probability 1 sets. Hence density of an attracting 
set in phase space does not mean that the corresponding attractor has dimension equal to that of the 
phase space: it could be substantially lower, see [GC95]. 
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distributions p. a obtained by letting the parameters a vary. We call such a collection 
an ensemble" . 

For instance a could be the average energy U of the system, the average kinetic energy 
T, the volume V, the intensity E of the acting non conservative forces, etc 

Non equilibrium thermodynamics can be defined as the set of relations that the vari- 
ations of the parameters a and of other average quantities are constrained to obey as 
some of them arc^ varied. In equilibrium the heat theorem is an example of such relations. 
In reversible non equilibria the fluctuation theorem (4.6) is an example. 

In non equilibrium systems the equations of motion play a much more prominent role 
than in equilibrium: in fact one of the main properties of equilibrium statistical mechanics 
is that dynamics enters only marginally in the definition of the statistical distributions 
of the equilibrium states. 

The necessity of a reversibility assumption in the fluctuation theorem already hints at 
the usefulness of considering the equations of motion themselves as "parameters" for 
the ensembles describing non equilibrium stationary states: we are used to irreversible 
equations in describing non equilibrium phenomena (like the heat equation, the Navier 
Stokes equation, etc) and unless we are able to connect our experiments with reversible 
dynamical models we shall be unable to make use of the fluctuation theorem. 

Furthermore it is quite clear that once a system is not in equilibrium and thermostat- 
ting forces act on it, the exact nature of such forces might be irrelevant within large 
equivalence classes: i.e. it might be irrelevant which particular "cooling device" we use 
to take heat out of the system. Hence one would like to have a frame into which to set 
up a more precise analysis of such arbitrariness. Therefore we shall set 

Definition 1: A stationary ensemble £ for a system of particles or for a fluid is the 
collection of SRB distributions, for given equations of motion, obtained by varying the 
parameters entering into the equations. 

It can happen that for the same system one can imagine different models. In this case 
we would like that the models give the same results, i.e. the same averages to the same 
observables, at least in some relevant limit. Like in the limit of infinite size in which 
the number N, the volume V and the energy U tend to infinity but N/V and U/V stay 
constant. Or in the limit in which the Reynolds number R tends to infinity in the case 
of fluids. 

This gives the possibility of giving a precise meaning to the equivalence of different 
thermostatting mechanisms. We shall declare 

Definition 2, (equivalence of ensembles): Two thermostatting mechanisms are equiv- 
alent "in the therm,odyna,m,ic limit" if one can establish a one to one correspondence 
between the elements of the ensembles £ and £' of SRB distributions associated with the 
two models in such a way that the same observables, in a certain class C of observables, 
home the same averages in corresponding distributions, at least when some of the param- 
eters of the system are sent to suitable limiting values to which we assign the generic 
name of "thermodynamic limit". 

In the following sections we illustrate possible applications of this concept. 
§6. Drude— Lorentz' electric conduction models. 

Understanding of electric conduction is in a very unsatisfactory state. It is usually 
based on linear response theory and very seldom a fundamental approach is attempted. 
Of course this is so for a good reason, because a fundamental approach would require 
imposing an electric field E on the system and, at the same time, a thermostatting force 
to keep the system from blowing up and to let it approach a steady state with a current 
Je flowing in it, and then taking the ratio Je/E (with or without taking also the limit 
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asE^ 0). 

However, as repeatedly mentioned, it is an open problem to study steady states out of 
equilibrium. Hence most theories have recourse to linear response where the problem of 

studying stationary non equilibria docs not even arise. 

The reason why this is unsatisfactory is that as long as we are in principle unable to 
study stationary non equilibria we are also in principle unable to estimate the size of the 

approximation and errors of linear response. 

In spite of many attempts the old theory of Drude, see [Be59], [Se87], seems to be among 
the few conduction theories which try to establish a conductivity theory based on the 
study of electric current at non zero fields. 

We imagine a set of obstacles distributed randomly or periodically and among them 
conduction electrons move, roughly with density of one per obstacle. 

The (screened) interactions between the electrons are, at a first approximation, ignored. 
The collisions between electrons and obstacles ("nuclei") will take place in the average 
after the electrons have traveled a distance A = {pa^)~^ if p is the nuclei density and a 
is their radius. 

Between collisions the electrons, with electric charge e, accelerate in the direction of an 
imposed field E incrementing, in that direction, velocity by 

Sv='-^ = ^^^ (6.1) 

where ks is Boltzmann's constant. At collision they are "thermalized" : an event that is 
modeled by giving them a new velocity of size v — ^ ksT /m and a random direction. 
The latter is the "thermostatting mechanism" which is a, somewhat rough, description 
of the energy transfer from electrons to lattice which physically corresponds to electrons 
losing energy in favor of lattice phonons, which in turn arc kept at constant temperature 
by some other thermostatting mechanism which prevents the wire melting. All things 
considered the total current that flows will be 

6^ d f 

^ pa?\/mkBT ^ ^ ^ 

obtaining Ohm's law. 

To the same conclusion we arrive by a different thermostat model. We imagine that 
the electrons move exchanging energy with lattice phonons but keeping their total energy 
constant and equal to N ks T: i.e. 12jLi mi^ — 'iNkBT/2, where fcs is Boltzmann's 
constant. There are several forces that can achieve this result 

we select the "Gaussian minimal constraint" force. ^ This is the force that is required to 
keep J^ rn-ij strictly constant and that is determined by "Gauss least effort" principle, 
see [Ga99a], ch. 9, appendix 4, for instance: as is well known this is, on the i-th particle, 
a force 

. dei eE-J2j ij . mE-N J . 

-ax,= ^-^^i - - -^M^ (6-3) 

If there are N particles and N is large it follows that J = N~^e^^ Xj is essentially 
constant, see [RuQQb], and each particle evolves, almost independently of the others, 
according to an equation: 



^ Not because it plays any fundamental role but because it has been studied by many authors and because 
it represents a mechanism very close to that proposed by Drude. We recall, for copleteness, that the 
effort of a constraint reaction on a motion on which the active force is / (with 3JV components) and a 
is the acceleration of the particles (with 3N components) and m is the mass is f ( a) = ( / —m a)^ /m; 
then Gauss' principle is that the effort is minimal if a is given the actual value of the acceleration, at 
fixed space positions and velocities. 
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between collisions, with a suitably fixed constant v. If wc imagine that the velocity of the 
particles between collisions changes only by a small quantity compared to the average 
velocity the "friction term" which in the average will be of order E"^ will be negligible 
except for the fact that its "only" effect will be of insuring that the total kinetic energy 
stays constant and the speeds of the particles are constantly renormalized. In other words 
this is the same as having continuously collisions between electrons and phonons even 
when there is no collision between electrons and obstacles. Hence the resulting current 
is the same (if N is large) as in (6.2). 

§7. Ensemble equivalence: the example of electric conduction theories. 

We have derived three models for the conduction problem, namely 

(1) the classical model of Drude, [Se87], in which at every collision the electron velocity 
is reset to the average velocity at the given temperature, with a random direction, c./.r. 
(6.1) and (6.2). 

(2) the Gaussian model in which the total kinetic energy is kept constant by a thermostat 
force 

fnXf=E— ~ i j + "collisional forces" (7.1) 

where SNksT is the total kinetic energy (a constant of motion in this model). The model 
has been widely studied and it was introduced by Hoover and Evans (see for instance 
[HHP87] and [EM90]). 

(3) a "friction model" in which particles independently experience a constant friction 

mi j = E_ — V x_^ + "collisional forces" (7-2) 

where v is a. constant tuned so that the average kinetic energy is eNksT 12. This model 
was considered in the perspective of the conjectures of ensemble equivalence in [Ga95] , 
[Ga96b]. 

The first model is a "stochastic model" while the second and third are deterministic: the 
third is "irreversible" while the second is reversible because the involution I{x^, v J = 
{x^, — Vi) anticommutes with the time evolution flow St defined by the equation (7.1): 

ISt = S-tl (as the "friction term" is odd under /). 

Let i^s,T be the SRB distribution for (7.1) for the stationary state that is reached starting 
from initial data with energy 3NkBT/2. The collection of the distributions ij,s,t as the 
kinetic energy T and the density S = N/V vary, define a "statistical ensemble" £ of 
stationary distributions associated with the equation (7.1). 

Likewise we call jls,^ the class of SRB distributions associated with (7.2) which forms 

an "ensemble" £. 

We establish a correspondence between distributions of the ensembles £ and £: we say 
that iJ,s,T and fls\i, are "corresponding elements" if 

6 = 5', T = j ^C^mx])jl5,u{dxdx) (7.3) 
Then the following conjecture was proposed in [Ga96b]. 

Conjecture 1: (equivalence conjecture) Let F be a "local observable", i.e. an observ- 
able depending solely on the microscopic state of the electrons whose positions is inside 
a fixed box Vo . Then, if £ denotes the local smooth observables 
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iV— >(x),JV/y=5 JV— »oo,Ar/y=5 

ifT and u are related by (7.3). 

This conjecture has been discussed in [Ga95], sec. 5, and [Ga96a], see sec. 2 and 5: and 
in [Ru99b] arguments in favor of it have been developed. 

Clearly the conjecture is very similar to the equivalence in equilibrium between canonical 
and microcanonical ensembles: here the friction v plays the role of the canonical inverse 
temperature and the kinetic energy that of the microcanonical energy. 

It is remarkable that the above equivalence suggests equivalence between a "reversible 
statistical ensemble", i.e. the collection £ of the SRB distributions associated with (7.1) 
and a "irreversible statistical ensemble", i.e. the collection £ of SRB distributions asso- 
ciated with (7.2). 

Furthermore it is natural to consider also the; collection £' of stationary distributions 
for the original stochastic model (1) of Drude, whose elements /i'^ ^ can be parameterized 
by the quantities T, temperature (such that | J2j n^i^ = ^NksT, and N/V = 6). This 
is an ensemble £' whose elements can be put into one to one correspondence with the 
elements of, say, the ensemble £ associated with model (2), i.e. with (7.1): an element 
/Lt'^ y G £' corresponds to fis^i, G £ HT verifies (7.3). Then 

Conjecture 2: If 1^5 € £ o,nd ^'g^ & £' are corresponding elements (i.e. (7.3) holds) 
then 

N^oo,N/V=5 N^oo,N/V=5 

for all local observables F G £. 

Hence we see that there can be statistical cqiiivalcncc between a viscous irreversible 
dissipation model and either a stochastic dissipation model or a reversible dissipation 
model, at least as far as the averages of special observables are concerned. 

The argument in [Ru99b] in favor of conjecture 1 is that the coefficient a in (6.3) is 
essentially the average J of the current over the whole box containing the system of 
particles, J = e J2j bence J should be constant with probability 1, at least if 
the stationary SRB distributions can be reasonably supposed to have some property of 
ergodicity with respect to space translations. 

§8. Entropy driven intermittency in reversible dissipation. 

A further argument for the equivalence conjectures in the above electric conduction 
models can be related to the fluctuation theorem: the quantity a{x) is also proportional 
to the phase space contraction rate a{x) = {3N — l)a{x). Therefore, denoting in general 
with a subscript + the SRB average (or the time average) of an observable, the probability 
that a{x) deviates from its average a+ = {3N — 1) a+ can be studied as follows. 

If the number N of particles is large the time scale tq over which a{Six) evolves will 
be large compared to the microscopic evolution rates, because (Ti(x) is the sum of the 
~ 6N rates of expansion and contraction of the ~ 6N phase space directions out of x 
(sometimes called the ^Hocal Lyapunov exponents").'^ 



The exact number of exponents depends on how many constants of motion the system has: for instance 
in the case of the conduction model (1) in §6 above the number of exponents is 6iV — 1 because the 
kinetic energy is conserved and the system has no other (obvious) first integrals. Furthermore one of such 
exponents is since every dynamical system in continuous time has one zero exponent (corresponding 
to the direction i of the flow). 
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Consider a large inimber m of time intervals Ii, L2. ■ ■ ■ , Im of size tq and let aj be the 
(average) value of a{Stx) for t & Ij. Then the fraction of the j's such that aj — C7+ ~ a+p 
will be proportional to 

7r,„(p):^e-°«(f) (8.1) 

and C{p) < C(l) if p 7^ 1. Since we can expect thatC(p) is proportional to N wc sec that 
the fraction of time intervals Ij in which aj ^ a+ will be exponentially small with N . 
For instance the fraction of time intervals in which cjj ~ — cr+ will be, by the fluctuation 
theorem 

g-(3JV-l)a+ro (82) 

In order that the above argument holds it is essential that A'' is large to the point that 
we can think that the time scale tq over which a{Stx) varies is much larger than the 
microscopic scales: so that we can regard tq large enough for the fluctuation theorem to 
apply. In this respect this is not really different from the previously quoted argument in 
[Ru99b] . However the change of perspective gives further information. 

In fact we get the following picture: A'' is large and for most of the time the (stationary) 

evolution uneventfully proceeds as if a{Stx) = (thus justifying conjecture 1). Very 
rarely, however, it proceeds as if a{Stx) ^ (7_|_, for instance with a{Stx) = — (7_|_: such 
"bursts of anomalous behavior'^ occur very rarely. But when they occur "everything else 
goes the wrong way" because, as discussed in detail in [Ga99c], while the phase space 
contraction is opposite to what it "should be" (in the average) then it also happens that 
all observables evolve following paths that are the time reversal of the expected paths, This 
is the content, see [Ga99c], of the following theorem which is quite close (particularly 
if one examines its derivation) to the Machlup-Onsager theory of fluctuation patterns 
(note that, however, it does not require closeness to equilibrium) 

Theorem (conditional reversibility theorem): If F is an observable with even (or odd), 
for sim,plicity, time reversal parity and if r is large then the evolution or "fluctuation 
pattern" (p{t) and its time reversal I(p{t) = ip{—t), t € [— ro/2, to/2], will be followed with 
equal likelihood if the first is conditioned to an entropy creation rate p and the second to 
the opposite —p. 

In other words systems with reversible dynamics can be equivalent to systems with 
irreversible dynamics but they show "intermittent behavior'^ with intermittency lapses 
that become extremely rare very quickly as ^ oo. Sometimes they can be really 
dramatic, as in the cases in which a ~ —cr+- alas they are unobservable just for this 
reason and one can wonder (see §9 below) whether this is really of any interest. 

§9. Local fluctuations and observable intermittency. 

As a final comment upon the analysis of the equivalence of ensembles attempted above 
we consider a very large system with volume V and a small subsystem of volume Vq which 
is large but not yet really macroscopic so that the number of particles in Vq is not too 
large, a nobler way to express the same notion is to say that we consider a "mesocopic" 

subsystem of our macroscopic system. 

Here it is quite important to specify the system because wc want to make use of as- 
pects of the equivalence conjectures that are model dependent. Therefore we consider 
the conduction models (2) or (3) of §5: these are models in which dissipation occurs 
"homogeneously" throughout the system. In this case we can imagine to look at the part 
of the system in the box Vq: if ji, . . . ,jNo are the particles which at a certain instant are 
inside Vq and Xj = f .(x ) are the equations of motion we can define 
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^Vo {x) = S^,, /i, {x) (9.1) 

i=l 

which is (by definition) the part of phase space contraction due to the particles in Vq- 
Since the part of the system inside the microscopically large but macroscopically small 
V{) can be regarded as a new dynamical system whose properties should not be different 
from the ones of the full system enclosed in the full volume V we may expect that the 
subsystem inside Vq is in a stationary state and the quantity has the same fluctuation 
properties as av, i-e. 

(1) {(^Vo)+ = Voa+, {av)+ = Va+ 

(2) ip) = e«f) - ^0 , w^ip) = ef(p) - ^ ^ ■ ^ 

where C)^+ are the same for V,Vo and p = fl!^/2^Vo{'^tx)/{(TvQ) dt or respectively 

p = J^!^^2^v{Stx)/{av) dt. Here avo is naively defined as the contribution to a 
coming from the particles in Vq. 

In other words in large stationary systems with homogeneous reversible dissipation phase 
space contractions fluctuate in an extensive way, i.e. they are regulated by the same 
deviation function C(p) (volume independent). 

This is very similar to the well known property of equilibrium density fluctuations in a 
gas of density p: ii V D Vq are a very large volume F in a yet larger container and Vq 
is a small but microscopically large {i.e. mesoscopic) volume Vq then the total numbers 
of particles in V and Vq will be N and Nq and the average numbers will be pV and pVo 
respectively. Then setting 

p={N-pV)/pV, or, p={No-pVo)/pVo (9.3) 
the probability that the variable p has a given value will be proportional to 

TT^{p) = Sp)v, 7r^«(p) = e«(f)^« (9.4) 
again with the same function C{p). 

This means that wc can observe (^{p) by performing fluctuaMons experiments in small 
boxes, ideally carved out of the large container, where the density fluctuations are not 
too rare. A "local fluctuation law" should hold more generally in cases of models in 
which dissipation occurs homogeneously across the system, like the above considered 
conduction models. 

The intuitive picture for the above "local fluctuation relation" inspired (and was sub- 
stantiated) a mathematical model in which a local fluctuation relation can be proved as 
a theorem: it has een discussed in [Ga99c], see also below. 

Going back to the conduction model we see that the intermittency phenomena discussed 
above can be actually observed by looking at the fluctuations of the contribution to phase 
space contraction due to a small subsystem. 

And such entropy driven" intermittency will be model independent for models which are 
equivalent in the sense of the previous sections provided the models used are equivalent 
and one of them is reversible. 

An extreme case is provided by models (1)%(3), §7, for electric conduction (conjectured 
to be equivalent, see §7). In fact at first the model (3), the viscous thermostat, might 
look uninteresting as, obviously, in this case 

a^{x) = 3Nu, a^°{x) = SNqu (9.5) 
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and jV has no fluctuations. 

However the equivalence conjecture makes a statement about expectation values of the 
same observable: hence we should consider the quantity &^°{x) = E_ ■ J_yg/J2j il'j 
wc should expect that its statistics with respect to an clement of the ensemble £' is 
the same as that of the same quantity with respect to the corresponding elements of the 
ensembles £, £. Hence in particular the functions C,{p) which control the large fluctuations 
of cF^{p) will verify 

~a-p) - ~ap)-p{i^'')+iv^ = Hp) - 3p^P = c{p) - ^^^p (9.6) 

where the first equality expresses the validity of a fluctuation theorem type of relation 
due to the fact that the small system, by the equivalence conjecture, should behave as a 
closed system; the second equality expresses a consequence of the equivalence conjecture 
between models (2) and (3) while the third is obtained by expressing the current via 
Drude's theory (again assuming the conjectures of equivalence 1,2 of §7). 

§10. Fluids. 

The chaotic hypothesis was originally formulated to understand developed turbulence, 
[Ru78] : it is therefore interesting to revisit fluid motions theory. 

The incompressible Navier Stokes equation for a velocity field u in a periodic container 
V of side L can be considered as an equation for the evolution in time of its Fourier 
coefficients u j. where the "mode" k has the form 27ri~^ n with n ^ and n an integer 
components vector.^ Furthermore Uk= k k ■ u k =0- If P is the pressure field 
and / a simple forcing we shall fix the ideas by considering / ( x ) = / e sin fc • x 
where is some prefixed mode and e is a unit vector orthogonal to fc j. 

The Navier Stokes equation is then 

u + u ■ d u ~ dp+ I + lyAu (10.1) 

and it is convenient to use dimensionless variables Uo,Po, ^ q, ^ ,t: so we define them as 
u{x,t)=fL^u~'^UQ{L~'^x,L~'^irt), |^=L~^^, T = L~'^ut 

p{x,t)=fLpo{L-^x,L-^ut), R'^^fU'v-^ (10.2) 

l{x,t) =f^^{L-^x) 

with max I (;9 pi = 1. The result, dropping the label and calhng again x,t the new 
variables ^ , t, is that the Navier Stokes equations become an equation for a divergenceless 
field u defined onV= [0, 1]^, with periodic boundary conditions and equations 

u+Ru-du=-dp+ip+Au, d-u=0 (10.3) 

with max | <^ | = 1. 

Equation (10.3) is our model of fiuid motion, where R plays the role of "forcing intensity" 
and the term Au represents the "thermostatting force". As R varies the stationary 
distributions /Zij which describe the SRB statistics of the motions (10.3) define a set £ 
of probability distributions which forms an "ensemble" . 

The mathematical theory of the Navier Stokes equations is far from being understood: 
however phenomenology establishes quite clearly a few key points. The main property 



The value n = is excluded because, having periodic boundary conditions, it is not restrictive to 
suppose that the space average of u vanishes (galilean invariance) . The convention for the Fourier 
transform that we use is 3i(2) = X^^ e*-— 
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is that if (10.3) is written as an equation for the Fourier components of u then one can 
assume that u^. = for | fc | > K{R), for some finite K[R). 

Therefore the equation (10.3) should be thought of as a "truncated equation" in mo- 
mentum space by identifying it with the equation obtained by projecting also u-du on 

the same function space. 

Slioiikl one develop anxiety about the mathematical aspects of the Navier Stokes equa- 
tion one should therefore think that an equally good model for a fluid is the mentioned 
truncation provided K{R} is chose large enough. 

The idea is that for K(R) = R'^, with k larger than a suitable kq the results of the 
theory, i.e. the statistical properties of become K-independent for R large. 
The simplest evaluation of kq gives kq = 9/4 as a consequence of the so called K41 

theory of homogeneous turbulence, see [LL71]. 

If (10.3) is a good model for a fluid when L is large then it provides us with an "ensemble" 
£ of SRB distributions (on the space of the velocity fields components u j. of dimension 
~ 8nK{Rf/3).^ 

We should expect, following the discussion of the statistical mechanics cases, that there 
can be other "ensembles" £ which are equivalent to £. 

Here R plays the role of the volume in non equilibrium statistical mechanics, so that R 
oo will play the role of the thermodynamic limit, a limit in which the effective number 
of degrees of freedom, ^ AttR^'^/S, becomes infinite. The role of the local observables 
will be plaid by the (smooth) functions F{u) of the velocity fields u which depend 
on u only via its Fourier components that have mode k with |fc| < B for some B: 
Fiu)=F{{uk}\k\<B)- 

We shall call £ the space of such observables: examples can be obtained by setting 
F{u) = I J e'-'- u{x)dx\'^ or F{u) = J tp{x)- u{x)dx where the function has only 
a finte number of harmonics, ip{x) = / ^ |<b e'-'- u (^) , etc. 

As in non equilibrium statistical mechanics we can expect that the equations of motion 
themselves become part of the definition of the ensembles. For instance one can imagine 
defining the ensemble £ of the SRB distributions jlv for the equations 

u+Ru-du=-dp+'tp_+v{u)i^u (10.4) 

called GNS equations in [Ga97a], or "gaussian Navier Stokes" equations, where ^{u) is 
so defined that 

E = [{duf dx/{2-Kf = 5^ fc' I « , 1^ (10.5) 

is exactly constant and equal to S. The equations (10.4) are interpreted as above with 
the same momentum cut off K{R) = R'^. 

An element /is of £ and one iir of SRB distributions for the two different dynamics 
(10.3) and (10.4), "correspond to each other" if 

S = j m{du){J^idufdx/{2nf)=^En (10.6) 

where /i/j e 5 is the SRB distribution at Reynolds number R for the previous viscous 
Navier Stokes equation, (10.3), and we naturally conjecture 

There are about 4TrK(R)^ /3 vectors with integer components inside a sphere of radius K(R), thus the 
number of complex Fourier components with mode label | fc | < K{R) would be 3 times as much, but 
the divergenceless condition leaves only 2 complex components for u j. along the two unit vectors 

orthogonal to k and the reality condition further divides by 2 the number of "free" components. 
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Conjecture 3 (equivalence GNS NS): If R ^ oo then for all local observables F G C it 
is ij,r{F) = /iHfl(-F) if (10.6) holds. 

It is easy to check that the GNS model "viscosity" v{u), having to be such that the 
quantity S in (10.5) is exactly constant must be 

Jy{^- Au — i?Au • [u ■ du)) dx 

and we see that while (10.3) is an irreversible equation the (10.4) is reversible, with time 
reversal symmetry given by 

Iu{x,t) = -u{x,t) (10.8) 

as one can check. 

More generally one may wish to leave the "Kolmogorov parameter" k as a free pa- 
rameter: in this case the SRB distributions will form an ensemble whose elements can 
be parameterized by R, k and the equivalence conjecture can be extended to this case 
yielding equivalence between ^ib.,k and /tH,re- This is of interest, particularly if one has 
numerical experiments in mind. 

If K > Ko then the value of k should be irrelevant: but if k < kq the phenomenology 
will be different from the one of the Navier Stokes equation and equivalence might still 
hold but one cannot expect either equation to have the properties that we expect for 
the usual Navier Stokes equations {i.e. in this situation one would have to be careful in 
making statements based on common experience). 

If we take k to be exactly equal to the value kq = 9/4 (i.e. if we take the ultraviolet 
cut-off to be such that, according to the K41 theory, for larger values it is needlessly large 
and for lower values it is incorrectly low and shows a phenomenology which will depend 
on its actual value) then we may speculate that the "attracting set" is the full phase 
space (available compatibly with the constraint S = S^). Therefore the divergence of 
the equations of motion, which is given by a rather involved expression in which only the 
first term seems to dominate at large R, namely 

a{u) = { k^)iy{u)-{ f A^-Audx){ f [{Auf- 

\k\<K{R) 

- RAu ■ (A(u -du)) - R{Au) ■ {Au) ■ (du) - RAu ■ {Adu)u+ (10.9) 



iy{u)Au ■ A'^ujdx) / / {Aufdx 
Jv 



will veriiy the fluctuation theorem, i.e. the rate function (^{p) for the average phase space 
contraction p = ^^^t ^ )dt/Ta+ will be such that C(— P) = C(p) ~ P<^+- 

If the chaotic hypothesis is valid together with the equivalence conjecture the validity 
of the fluctuation relation can be taken as a criterion for determining /t: it would be 
the last K before which the fluctuation relation between C,{p) and (,{—p) holds. However 
this conclusion can only be drawn if the attracting set in phase space is the full ellipsoid 
S = Sij at least for K{R) = i^"". 

The latter property might not be realized: and in such case the fluctuation theorem 
does not apply directly, although the equivalence conjectures still hold. In fact one can 
try to extend the fluctuation theorem to cover reversible cases in which the attracting set 
is smaller than the full phase space left available by the constraints. In such cases under 
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suitable geometric assumptions, [BG97] and the earlier work [BGG97], one can derive a 
relation like 

((-p) = ({p) -pa+'d, < < 1 (10.10) 

where i? is a coefficient that can be related to the Lyapunov spectrum of the system, 
c./.r. [BG97], [Ga97a]. In fact numerical work to check the theory proposed in [Ga97a] is 
currently being performed (private communication by Rondoni and Segre) with not too 
promising results which, optimistically, can be attributed to the fact that the ultraviolet 
cut off is too small due to numerical litmitations: clearly there is more work to do here. 
The preliminary numerical results give, so far, the somewhat surprising linearity in p but 
with a slope that, although of the correct order of magnitude, seems to have a value that 
does not match the theory within the error bounds. 

Coming back to the Navier Stokes equation we mention that we may imagine to write 
it as (10.3) hut with the different constraint 

U= I u^rfs = const (10.11) 
Jv 

rather than (10.5). 

This case has been considered in [RS99] and the multiplier v{u) is in this case 

- - \k\<K{R) 

and we can (almost) repeat the above considerations and equivalence conjectures. This 
constraint is a gaussian constraint that U is constant obtained by imposing its constancy 
on the Euler evolution via Gauss' principle with a suitable definition of the notion of 
"constraint effort" (this notion is not unique, see [Ga97a] for another definition) and we 
do not discuss it here to avoid overlapping with §12 below. 

Thd intuitive motivation for the equivalence conjectures is that for large R the phase 
space contraction (j{u) and the coefRcient i'{u)'^ are "global quantities" and depend 
on the global properties of the system {e.g. (t(m) is the sum of all the local Lyapunov 
exponents of the system whose number is 0{K{R)^)): they will "therefore" vary over 
time more slowly than any time scale of the system and can be considered constant. 

The argument is not very convincing in the case of the equations with the constraint 
(10.11) because the cr( m ) in (10.12) is proportional to Jy ^fi ■ u dx which clearly depends 
only on harmonics of u with k small, i.e. it is a "local observable". Note that this does 
not apply to the GNS equations with the constrained vorticity S, (10.6) where the "main" 
contribution to (t( w), see (10.7), comes from the term proportional to R which contains 
all harmonics. Therefore the result in [RS99] about the equivalence between the GNS 
equations, (10.4) with the constraint (10.5), and the equations with constraint (10.11) is 
interesting and puzzling: it might be an artifact of the smallness of the cut off that one 
has to impose in order to have numerically feasible simulations. 

Finally cr+(li)/cr+, i.e. essentially i'{u)/i'+ will fluctuate taking values sensibly differ- 
ent from their average value 1, at very rare intervals of time: but when such fluctuations 
will occurr one shall see "bursts" of anomalous behavior: i.e. the motion will be ^Hnter- 
mittenf as in the case discussed in non equilibrium statistical mechanics. 



Which in the case (10.9) are simply proportional and in the case of (10.4) they are related in a more 
involved way, see (10.8), (10.9), but which are still probably proportional to leading order as K — > oo. 
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11. Entropy creation rate cind entropy driven intermittency. 

Of course if R is large the number of degrees of freedom is large and intermittency on 
the scale of the fluid container will not be observable due to its extreme unlikelyhood 
(expected and quantitatively predicted by the fluctuation theorem). 

Therefore we look also here, in fluid motions, for a local fluctuation relation. Fluids seem 
particularly suitable for verifying such local fluctuations relations because dissipation 
occurs homogeneously, i.e. friction strength is translation invariant. 

This implies that we can regard a very small volume Vq of the fluid as a system in itself 
(as always done in the derivation of the basic fluid equations, e.g. see [Ga97b]) and we 
can expect that the phase space contraction due to such volume elements is simply cr( w), 
given by (10.9) or (10.12) ( "equivalently" because of our equivalence conjectures) with 
the integrals in the numerator and denominator being extended to the volume Vb rather 
than to the whole box, and expressing (essentially by definition) the "local phase space 
contraction" crvo(w). 

Then p = J^^^^^.^avu{Siu)/{(TVa)+ will have a rate hmction C(p) which will verify, 
under the same assumptions as in (10.10), a large deviation relation as 

a-P) = ap)-p{<Tvo)+^ (11.1) 

for some d: as mentioned the theoretical value of this slope seems currently inaccessible 
to theory (as the theory proposed in [BG67] , [Ga97a] may need substantial modifications, 
c.f.r. comment following (10.10)). The (crvo) and C(p) will be proportional to Vq: ({p) = 
Vo C{p) with a Vo-independent C{p). Note that C(p) depends also on R. 

The small volume element of the fluid will therefore be subject to rather frequent varia- 
tions: in spite of ({p) being proportional to Vq, because now Vq is not large. The conse- 
quent intermittency phenomena can therefore be observed. And as in ^9 once the phase 
space contraction is intermittent all properties of the system show the same behavior. 

And in fact intermittency in observations averaged ove a time span r will appear with a 
time frequency of the form e^" (C(p)-C(i)) t^j^g quantity p can be interpreted as a measure 
of the ''strength of intermittencij' observable in easurements averaged over a time r 
because as noted in §9 and in [Ga99b] the size of p controls the statistical properties 
of "most" other observables. Therefore the function Q{p) (hence C,{p) ) might be directly 
measurable and it should be rather directly related to the quantities that one actually 
observes in intermittency experiments. And the difference Q{p) — ({—p) can be tested for 
linearity in p as predicted by the analysis above. 

Note that in an extended system the volume V is much larger than Vq and we shall 

see "for sure" intermittency (for observables averaged over a time r) of strength p in a 
region of volume Vq somewhere within a volume W such that 

At this point it seems relevant to recall that it is rather heatedly being debated whether 
the name of ''entropy creation rate" that some authors (including the present one) give 
to the phase space contraction rate is justified or not, see [An82]. The above properties 
not only propose the physical meaning of the quantity p and bring up the possibility 
of measuring its rate function ^(p) in actual experiments but also provide a further 
justification of the name given to a as "entropy creation rate" and fuel for the debate 
that inevitably the word entropy generates at each and every occurrence. 



(11.2) 
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§12. Benard convection, intermittency and the Ciliberto— Laroche experi- 
ment. 



A very interesting attempt at checking some of the above ideas has been made recently 
by Cilberto and Laroche in an experiment on real fluids which has been performed with 
the aim of testing the relation (11.1) locally in a small volume element, [CL98]. By "real" 
we mean here non numerical: a distinction that, however, has faded away together with 
the XX th century but that some still cherish: the system is physically macroscopic 
(water in a container of a size of the order of a liter). 

This being a real experiment one has to stretch quite a bit the very primitive theory 
developed so far in order to interpret it and one has to add to the chaotic hypothesis 
other assumptions that have been discussed in [BG97], [Ga97a] in order to obtain the 
fluctuation relation (10.12) and its local couterpart (11.1). 

The experiment attempts at measuring a quantity that is eventually interpreted as the 
difference ^(p) — C{~P): by observing the fluctuations of the product -du^ where is the 
deviation of the temperature from the average temperature in a small volume element A 
of water at a fixed position in a Couctte flow and is the velocity in the z direction of 
the water in the same volume element. 

The result of the experiment is in a way quite unexpected: it is found that the function 
C(p) is rather irregular and lacking symmetry around p — 1: nevertheless the function 
C(p) — C(~P) seems to be strikingly linear. As discussed in [Ga97a], predicting the slope of 
the entropy creation rate would be difficult but if the equivalence conjecture considered 
above and discussed more in detail in [Ga97a] is correct then we should expect linearity 



In the experiment of [CL98] the quantity {iu^ did not appear to be the divergence of 
the phase space volume simply because there was no model proposed for a theory of the 
experiment. Nevertheless Ciliberto-Laroche select the quantity j^-du^ dx on the basis 
of considerations on entropy and dissipation so that there is hope that in a model of the 
flow this quantity can be related to the entropy creation rate discussed in §10, §11. 

Here we propose that a model for the fluid, that can be reasonably used, is Rayleigh's 
model of convection, [Lo63], [LL71] and [Ga97b] sec. 5. An attempt for a theory of the 
experiment could be the following. 

One supposes that the equations of motion of the system in the whole container (of 
linear size of the order of 30 cm) are written for the quantities x, z, ??, u in terms of the 
height H of the container (assumed to be a horizontal infinite layer), of the temperature 
difference between top and bottom 6T and in terms of the phenomenological "friction 
constants'^ X of viscosity, dynamical thermal conductivity and of the thermodynamic 
dilatation coefficient a. We suppose that the fluid is 3 dimensional but stratified, so 
that velocity and temperature fields do not depend on the coordinate y, and gravity is 
directed along the z-axis: g = g e, e — (0,0,-1). The temperature deviation i? is 
defined as the difference betwen the temperature T{x, y, z) and the temperature that the 
fluid would have at height z in absence of convection, i.e. Tq — z 5T / H HTq is the bottom 
temperature. 

In such conditions the equations, including the boundary conditions (of fixed temper- 
ature at top and bottom and zero normal velocity at top and bottom), the convection 
equations in the Rayleigh model, see [Lo63] eq. (17), (18) where they are called the 
Saltzman equations, and [Ga97b] §1.5, become 



ofC(p)-C(-p)- 



d ■ u = 




12.1 



u + u ■ d u = vAu — a'&g — dp', 



(12.1) 



m 
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12.3 



12.4 



12.5 



12.6 



The function p' is related to the pressure p: within the approximations it is p = Pa — 
Pogz + p'. We shall impose for simplicity horizontal periodic boundary conditions in x, y 
so that the fluid can be considered in a finite container V of side a for some a > prefixed 
(which in thcoriginal variable would correspond to a container of horizntal size aH). 
It is useful to define the following adimensional quantities 



ad 

12.2 a 



= u' = iV^H^r' u (12.2) 



2 gH^aST 

K = 5 , Hpr = — 

and one checks that the Rayleigh equations take the form 
u + Ru ■ d u = Au — i?i?e — dp, 



d-u=0 (12.3) 

u,(0) =M,(1) = 0, ^(0) = i?(l) = 0, 



/ Uxdx = I Uydx = 
Jv Jv 



where we again calH, x,y,z, u, i? the adimensional coordinates r,^, 77, C, u°, in (12.2). 
The numbers R,Rpr are respectively called the Reynolds and Prandtl numbers of the 
problem: Rpr ='~ 6.7 for water while i? is a parameter that we can adjust, to some 
extent, from up to a rather large value. 

According to the principle of equivalence stated in [Ga97a] here one could impose the 
constraints 

[ (u^ + -^d^)dx =C (12.4) 

Jv ^Pr 

on the ^^frictionless equations" {i.e. the ones without the terms with the laplacians) and 
determine the necessary forces via Gauss' principle of minimal effort, see footnote ^ 
and [Ga96a], [Ga97a]. We use as ejfort functional of an acceleration field a and of a 
temperature variation field s the quantity 

^(a,s)'= ((a + dp- l),{-A)-\a + dp- /))+ (12.5) 
+ ((s-<^),(-A)-i(s- V)) with 

and require it to be minimal over the variations S{x) of a = ^ and t{x) ofs = ^ 
with the constraints that for all x it is 9 • ^ = 0, besides those due to the boundary 
conditions. The result is 

d-u=0 

u + Ru ■ du = R-d e — dp' -\- Tth 

'd + Ru ■ d-d = Ru' + Xth (12-6) 



= '&{H), / Uxdx = / Uydx = 
Jv Jv 
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12.7 



where the frictionless equations are modified by the therm,ostats forces Tiih-^ih'- the 
latter impose the nonholonomic constraint in (12.4) with the effort functional defined by 
(12.5). Looking only at the bulk terms we see that the equations obtained by imposing 
the conlstraints via Gauss' principle become the (12.3) with coefficients in front of the 
Laplace operators equal to VG^^GRpr, respectively, with the "gaussian multipher" 
being an odd functions of u, see [Ga97a]: setting Cv{lk,'>^) — Jy {{du)^+Rpl.{di))'^)dx 
one finds 

UG = Cv{u,'d)-'^R{l+Rpl) I u'^dx (12.7) 

Jv 

And the equations become, finally 
d-u=0 

u + Ru ■ d u = R'd e — dp' + vq 



12.8 '& + Ru-d'd = Ru^ + UG^^'& 



(12.8) 



^(0) = = '(?(i?), j u^dx = j Uydx = 



If one wants the equivalence between the ensembles of SRB distributions for the equation 
(12.8) and for (12.3) one has to tune, [Ga97a], the value of the constant C in (12.4) so 
that the time average value {i'g)+ of vg is precisely the physical one: namely {vg) = 1 
by (12.3). This is (again) the same, in spirit, as fixing the temperature in the canonical 
ensemble so that it agrees with the microcanonical temperature thus implying that the 
two ensembles give the same averages to the local observables. 

The equations (12.8) are time reversible (unlike the (12.3)) under the time reversal map: 



{u,i)) = {-u,d) (12.9) 

and they should be supposed, by the arguments in [Ga97a] and §10,11: "equivalent" to 
the irreversible ones (12.3), 

The (12.8) should therefore have a "divergence" (j{u, i9) whose fluctuation function ({p) 
verifies a linear fluctuation relation, i.e. (^{p) — C{~P) should be linear in p. Note that 
the divergence of the above equations is proportional to ug if one supposes that the 
high momenta modes with | fc| > K(R) = R'^ with k suitable can be set equal to so 
that the equation (12.8) becomes a system of finite differential equations for the Fourier 
components of u,'&. 

For instance the Lorenz' equations, [Lo63] see also §17 of [Ga97b], reduced the number of 
Fourier components necessary to describe (12.3) to just three components, thus turning 
it into a system of three differential equations. 

Proceeding in this way the divergence of the equations of motion can be computed as a 
sum of two integrals one of which proportional to i^g in (12.7). If instead of integrating 
over the whole sample we integrate over a small region A, like in the experiment of 
[CL98], we can expect to see a fluctuation relation for the entropy creation rate if the 
fluctuation theorem holds locally, i.e. for the entropy creation in a small region. 

As for the cases in §11 this is certainly not implied by the proof in [GC95]: however 
when the dissipation is homogeneous through the system, as it is the case in the Rayleigh 
model there is hope that the fluctuation relation holds locally because "a small subsystem 
should be equivalent to a large one". As noted in §9 the actual possibility of a local 
fluctuation theorem in systems with homogeneous dissipation has been shown in [Ga99c] , 
after having been found through numerical simulations in [GP99] , and this example was 
relevant because it gave us some justification to imagine that it might apply to the present 
situation as well. 
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The entropy creation is due to the term R j^u^'O dx_/C/^{u ,%')), where A is the region 
where the measurements of [CL98] are performed, hence we have a proposal for the 
explanation of the remarkable experimental result. Unfortunately in the experiment 
[CL98] the contributions not explicitly proportional to R to the entropy creation rates 
have not been measured nor has been the C in (12.7) which also fluctuates (or might 
fluctuate). In any event they might be measurable by improving the same apparatus, so 
that one can check whether the above attempt to an explanation of the experiment is 
correct, or try to find out more about the theory in case it is not right. If correct the 
above "theory" the experiment in [CL98] would be quite important for the status of the 
chaotic hypothesis. 

§13. Conclusions. 

The chaotic hypothesis promises a point of view on non equilibrium that has proved so 

far of some interest. Here wc have exposed the basic ideas and attempted at drawing 
some consequences: admittedly the most interesting rely on rather phenomenological and 
heuristic grounds. They are summarized below. 

(1) The definition of nonequilibrium ensembles with the proposal that out of equilibrium 
also the equation of motion should be considered as part of the definition of ensemble. 
This is take into account that while in equilibrium the system is uniquely defined by 
its microscopic forces and constituents in non equilibrium it is not so. Systems must be 
put in contact with thermostats if we want them to become stationary after a transient 
time. And (for large systems) there may be several equivalent ways of taking heat out 
of a system, i. e. several thermostats, without affecting the properties of stationary state 
that is eventually reached by the system itself. 

(2) Equivalence of ensembles has the most striking aspect that systems which evolve with 
equations that are very different may exhibit the same statistical properties. In particular 
reversible evolutions might be equivalent to non reversible ones, thxis making it possible 
to apply results that require reversibility, in particular the fluctuation relations, to cases 
in which it is not valid. 

(3) An interpretation of the quantity p that intervenes in the fluctuation theorems in 
terms of an intermittency phenomenon and as a further quantitative measure of it. 

(4) The possibility of applying the theory to strongly turbulent motions was the origin of 
the Ruellc's principle that evolved into the chaotic hypothesis: therefore not surprisingly 
the ideas can be applied to fluid dynamics. We have discussed a possible approach. The 
approach leads again to a proposal for the theory of certain intermittency phenomena 
which appear quantitatively related to entropy creation fluctuations. 

(5) The possibility of measurement of the rate C(p) leads to a a possible prediction of the 
spatial frequency of internittent events of strength p or, as I prefer, with entropy creation 
ratep (sec (4) above, (11.2) and §12). This seems testable in concrete experiments (both 
real and numerical). 

(6) We have used the results in (2)% (8) to hint at an interpretation of the experiment 
by Ciliberto and Laroche on Benard convection in water. 

Although the theory is still at its beginning and it migh turn out to be not really of 
interest it seems that at this moment it is worth trying to test it both in its safest, c.f.r. 

§2%§8, and in its most daring, c.f.r. §9%§12, predictions. 
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